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Abstract 

The dynamical degenerate four-wave mixing is studied analytically in detail. 
By removing the unessential freedom, we first characterize this system by a lower- 
dimensional closed subsystem of a deformed Maxwell-Bloch type, involving only 
three physical variables: the intensity pattern, the dynamical grating amplitude, 
the relative net gain. We then classify by the Painleve test all the cases when 
singlevalued solutions may exist, according to the two essential parameters of the 
system: the real relaxation time t, the complex response constant 7. In addition 
to the stationary case, the only two integrable cases occur for a purely nonlocal 
response (3? (7) = 0), these are the complex unpumped Maxwell-Bloch system and 
another one, which is explicitly integrated with elliptic functions. For a generic 
response (^(7) 7^ 0), we display strong similarities with the cubic complex Ginzburg- 
Landau equation. 

Keywords: exact solutions, four- wave mixing, Painleve test, singularity analysis, solitary 
wave solutions. 
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1 Introduction 

The wave self-action by the degenerate mixing in a nonlinear medium involves three 
simultaneous processes: the interference of waves, the recording of the dynamical grating 
by an interference pattern, and the wave diffraction by the grating. This process is now 
the basic technique of important practical applications in real time holography, including 
optical phase conjugation, holographic interferometry, novelty filters, all-optical signal 
processing, etc [22 1 [151 fTTj . 

During the wave mixing, the self-diffraction of waves is governed by a self-consistent 
set of five equations for five complex amplitudes Aj,j = 1, 2, 3, 4 and £ , see e.g. [22] 

d z Ai = -i£A 2 , d z A 2 = i£A u d z A 3 = -i£A 4 , 3 Z A 4 = i£A 3 , (1) 

d t £ = 7 I m - -, (2) 
r 

I m = A{A 2 + A 3 A 4 , (3) 

where ([T]) is the coupled wave system for slow variable amplitudes Aj(z,t) [2S], © is the 
evolution equation of the grating amplitude £ with a rhs including the grating gain and 
the grating relaxation, ([3]) is the relevant interference pattern of the interacting waves. In 
our notation, bar denotes complex conjugation, d partial derivation, r is a real constant. 
It must be emphasized that the response constant 

7 = \l\e i9 (4) 

is complex. We will use the terms "local" and "nonlocal" response to describe the phase 
shift between the index grating £ and the interference pattern I m . In the case of a purely 
nonlocal response (7 purely imaginary), an energy transfer occurs between the interacting 
waves, whereas a local response (7 real) is characterized by an exchange of the phases of 
the waves [22] . In particular, the complex value of the coupling coefficient £ is an essential 
feature for the existence of soliton-like solutions. 

Apart from t and r, all variables are assumed dimensionless, after normalizing the 
physical variables A'j, z', 
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where k$ is the amplitude of the wave- vector in the free space, Io is the total input intensity 
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'0 



£4 

3=1 



constant, Ij 



A', 



(6) 



We restrict here to the so-called degenerate four- wave mixing (the four frequencies are 
identical), in the transmission geometry and in two space dimensions, 



kj kj^x&x kj^ z c Z} j 1, 2, 3, 4, 
k\ — k 2 = £4 — fc 3 = K, 

(e x and e z are unit vectors, K is the grating vector). 
So far, there exist two main analytic results: 



(7) 
(8) 



for 7 purely imaginary (purely nonlocal response) and in the stationary regime, a sech 
profile grating amplitude [15] ; 

when the phases of each Aj are independent of z, a parametric representation of the 
five amplitudes also restricted to a purely imaginary 7 [TH HI [5] , 

£ = (d z u) e^ e , 7 = «7nl, 7nl real, 

At = f 12 sin( Sl2 ( M - C 12 ))e*V\ A 2 = / 12 cos(s 12 ( M - C 12 ))e i( P*, 

i<Pa 



3J(7) = : < 



A 4 = -/ 43 sin(s 43 (w + C 43 ))e^ 4 , A 3 = / 43 cos(s 43 (w + C 43 ))e 

71 

tft - <p 2 - (p e + - = n 12 TC, S12 = (-1)™ 12 , 

V?4 ~ ~ ¥e + - = n 43 7l, S 43 = (-1)" 43 , 



(9) 



±J(<Pe ~ 7T/2) / f 2 



f l2 sin 2(u - C12) - / 43 sin 2(tt + C 43 ; 



. ™12>™43 G 

in terms of the real solution u of a damped sine-Gordon equation [TSJ IU E], 



w 2 t H — m 2 — Ksin(2ii + a) = 0, i^e 2 



7nl 



T 



(AV 



-2iCi; 



2 „2iC 43 



/43 e 



(10) 



The representation displays the invariance (1, 2, 3, 4, <9 Z , it) — > (4, 3, 2, 1, — d z , —u) 
and depends on six arbitrary real functions of t (f\ 2 , / 43 , Ci 2 , C 43 and the values of 
(fi + ip 2 and y? 4 + y5 3 ) an d one arbitrary real constant (the phase <p e ). The stationary 
sech solution [18] is then represented by [5] (see Eq. ff23l below), 



tew 



D 2A;(2; - z ) 



(11) 



In the present article we classify all cases when the system admits solutions with a 
singlevalued dependence on the initial conditions, and, with one major exception, we 
integrate all these cases. This major exception, left for future work, presents analogous 
difficulties to the search, in the complex cubic Ginzburg-Landau equation (CGL3), 

iA t +pA xx + q\A\ 2 A-i 1 A = 0, Pqi ^0, (A,p, q) E C, 7 G K, (12) 

for source [3], pulse [23] or front [21] solutions. 
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2 The intrinsic four- wave mixing, a deformed Maxwell- 
Bloch system 

The ten-dimensional system ([I])-® is invariant under any time-dependent rotation in the 
space {Ai, A%, A±, A3} which preserves the interference pattern ([3]). In order to remove 
this five-parameter unessential freedom, let us apply repeatedly the derivation operator 
d z , starting from the interference pattern ([3]), until a closed system has been obtained. 
This process ends after two steps and results in the intrinsic system 

£ 

d z I m = -iSI d , dzh = -2i€I m + 2i£I m , d t £ = ~/I m , (13) 

r 

admitting the first integral 

4|J m | 2 + il = K(t), K arbitrary. (14) 

The real field 1^ which is thus naturally introduced, 

I d = -\A X \ 2 + \A 2 \ 2 -\A^ + \A A \\ (15) 

has a natural interpretation: this is the relative net gain, therefore the four-wave mixing is 
characterized by three intrinsic variables: the intensity pattern J m , the grating amplitude 
£ and the relative net gain 

In previous integration methods [2] for the four-wave mixing, one would mainly look 
for the wave amplitudes Aj from some nonlinear system. Thanks to the existence of 
the above intrinsic system, the integration, whether analytic or numerical, now becomes 
systematic and involves two steps, 

1. Integration of the nonlinear intrinsic system (lT3"j) ; 

2. Knowing the grating £, integration of the two-dimensional linear system 

d z X = -i£Y, d z Y = -i£X; (16) 

indeed, given two linearly independent solutions (X,Y) = (Xi, Yi), (X%, Y 2 ), the 
general solution of (pQ) is defined in matrix form by 




in which the eight integration constants aij,b{j : constrained by the three relations 
(J3j) and ( TL5l) . depend on five arbitrary parameters according to relation (16T1) in the 
Appendix. 

The above system ([TBI is very similar to another classical system of nonlinear optics, 
the pumped Maxwell-Bloch system, which is an integrable system defined in complex 
form as [8] 

d x p = Ne, d x p = Ne, d x N = -{pe + pe)/2 + 4s = 0, d T e = p, d T e = p, (18) 
with s a real constant (the system is "pumped" when s is nonzero). 
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In fact, there is only one situation when the intrinsic four- wave mixing system f|T3|) 
and the pumped Maxwell-Bloch system ( 1T8]) can be identified. This occurs when, at the 
same time, the four-wave mixing model is undamped (r = +00) and has a purely nonlocal 
response (9?( 7 ) = 0), while the Maxwell-Bloch system is unpumped (s = 0). After this 
identification, 

i=o, S(T) =o, s =o: ±=«= m± = ?m = m = =^ = ™, (19) 

t X 1 p p JM e e 

the undamped, purely nonlocal response four-wave mixing model admits all the solutions 
of the unpumped complex Maxwell-Bloch system. 

The undamped case (relaxation time r = +00) physically means the recording of 
a permanent grating. In optics that can be, for example, the permanent holographic 
memory realized in nonlinear media. 

For practical computations, it may be advisable to eliminate i m from the grating 
evolution (T5]) and to equivalently consider the three-dimensional fifth order closed system, 

' | 7 | 2 <9 2 i d - 2i 1 £{d t £ + E/t) + 2i^E{d t £ + S/r) = 0, 
< (d K d t + ~d g )£ + i'y£I d = Q, (20) 
k 4| 7 |- 2 |d t £ + £/r\ 2 + I 2 d = K(t), K arbitrary. 

The sequel will display the crucial role of the third intrinsic variable (the relative net 
gain Id) to perform the explicit analytic integration whenever it is possible. 



3 The stationary case: general solution 

When the amplitudes are independent of the time t, the integration can be performed 
completely. The intrinsic system (fTHl) - (fHj) for I m ,I d ,£ reduces to 



1 -^m " 

4\I m \ 2 + li = K 



i£I d , —J d = -4|7|r(sin#)|i m | 2 , £ = ^rl n 
az 



(21) 



in which the first integral K is independent of t, therefore id obeys a first order ordinary 
differential equation (ODE) of the Riccati type, 



— i d = |7|r(sins)(i d 2 -i0. 



(22) 



The case 7 real is uninteresting for it involves no energy exchange and the intensities 



\£\ 



, |i m | , id are all constant. 

For 7 nonreal, the nonlinear intrinsic system ([2T]) admits the general solution 



id 
K 



k tanh kz 
| 7 |r sin g ' 



£ = ^rl a 



2 sing 



(k sech kz 



l-icotgg 



(23) 



\j7|TSing / 

in which k, Zq, <^o are constants of integration, with z — z written for shortness as z. 
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These bright profiles for \£\ 2 and |/ m | 2 extrapolate the solution of Ref. [5] which was 
restricted to 7 purely imaginary. 

The amplitudes are found by noticing that each variable A,£ _1//2 , j = 1,4 and A,£ -1 / 2 , 
j = 2,3 obeys a second order linear ODE with constant coefficients. The result is 



7 i K : < 



A l = {ksedikz) {l - lcotg9),2 e +iv °- ial2 ( 

A 3 = (ksechkzf +lcotg9)/2 ^ 0+i9 / 2 
A4 = (k sech kz) 



a±2 c_ + bu s_ 



A 2 = {ksechkzy^ l ^ y,, 'e-^ 0+l ^{ -a 12 s + + b 12 c H 



(l-icotgff)/2 ^ + i ipo -i g /2 i 



I-B34S- +a 34 c_; 

c± = cosh(l ± i cotg (7) — , s± = sinh(l ± i cotg g) — , 



(24) 



in which the conditions that Aj be complex conjugate of A, requires the four complex 
constants a 12 , &12, ^34, -B34 to be represented as 



' ai2 = Rcos\e iai \ 612 = Rcosfie 1 ? 12 , A u = RsmXe~ ia3 \ B u 

d2 fc sin(a 34 - P&) 

2K =- — j , — -r = — tan A tan a. 

7 r sing sm(ai2 - /3 12 ) 



R sin fie 



(25) 



The five additional constants of integration, chosen to be all real, are A,/i, a i2 + Pn, 
«34 + P34, and for instance an — /3u + a 34 — ^34 . 



4 Determination of the cases of singlevaluedness 

In the nonstationary case, the only existing analytic result, valid for a purely nonlocal 
response (^(7) = 0) and recalled in the introduction, is the parametric representation 
of the five complex amplitudes in terms of the solution u of the damped sine-Gordon 
equation (fTOj) . Rather than looking for solutions of this damped sine-Gordon equation, 
which would only concern the case $£(7) = 0, let us investigate the question of finding 
singlevalued closed form solutions, by applying the Painleve test [13] in order to detect 
all obstacles to singlevaluedness. 



4.1 The Painleve test 



For the basic notation (singular manifold variable (p, expansion variable x, auxiliary func- 
tions S,C), we refer to detailed lecture notes [TO] . 

Near a noncharacteristic (i.e. d z d t 7^ 0) movable singular manifold, as shown in our 
preliminary article [TTJ, the amplitudes have the leading order, 



A k ~ 
A k ~ 



a-kX 



-l+ib 
-1-ib 



A k ~ 
A k ~ 



£ 



%X 



-1+2* 



k 
k 



1,4, 
2,3, 



Jm ~ J m,oX , J d 



(26) 



,-2-2i& 



in which b is anyone of the two real constants defined in terms of 7 by 

(2b 2 — 1) cosg + 36 sin g = 0, = arg7. 



(27) 
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The leading coefficients depend on the nonzero auxiliary function C(z,t) and four arbi- 
trary complex functions A, /i,pi2,P43 of (z,t), 

ai = NX pi2 cosh /x, 6 2 = —NX p\ 2 x cosh /x, 

a 4 = NX P43 sinh /x, 63 = iVA PJ3 1 sinh /x, 

a 2 = A^A _1 pi2 cosh/i, 61= A^A -1 ^ 1 cosh/x, 
-1. 



a 3 = 

go = 

-^m,0 

N 2 -- 



—NX x p43 sinh /x, 64 = A 7 " A p 43 sinh /x, 
-i(l - i6)A 3 , r = + i6)A" 2 , 

= -iV 2 A 2 , TT 



c 

ItI 



'm,0 



A^A^, J,, 



-2AT 2 



1 - 26 2 )sinc; + 36cosc7) , C ^ 0. 



(2? 



The Fuchs indices of the linearized system only depend on the value of 6; for the 
ten-dimensional system ([I])-®, these are [TT] 



j = -1,0,0,0,0,2,2,2, 



5 ± Vl - 486 2 



(29) 



For the intrinsic five-dimensional system (1131) . the indices are 



5±y/T^48F 



(30) 



then the linear system (fl~6]) admits the Fuchs indices 

J = 0,2. 



(31) 



The diophantine condition that all Fuchs indices be integer therefore only admits the 
solution 6 = 0, 3^(7) = corresponding to a purely nonlocal response of the medium. 

In order to compute the necessary conditions for the absence of movable logarithms 
arising from the integer Fuchs indices, one can handle equivalently either the ten-dimensio- 
nal nonlinear system ([I])-©, or the five- dimensional nonlinear system (fT3"j) followed by 
the two-dimensional linear system (jl~6l) . One must distinguish 6 = from 6^0, and 
it is useless to test the quadruple index (because the leading order already introduces 
four arbitrary functions) and the index 4 (because of the existence of the singlevalued 
first integral K(t), Eq. (1141) ). In the generic situation 6 ^ no movable logarithm is 
detected at the triple index 2. In the nongeneric situation 6 = 0, for instance with the 
five-dimensional system ([TBI , two such necessary conditions Qj = are generated, at the 
Fuchs indices j ' = 2 and 3, 



9fc(7) = : 




r- 1 (C t + CC Z - (2/r)C) = 
1 A tt + (2/r)A t - 2CA 



T 



zt 



- C 2 A„ 



0, A 



(32) 



and no additional condition arises from the Fuchs index 2 of the linear system ([16]) . 

Remark. The analysis of the damped sine-Gordon equation (fTOl) only generates the 
condition Q 2 = [11] , since the condition Q 3 = which involves the phases of the complex 
amplitudes is then identically satisfied. 

A first solution to the conditions (I3"2"|) is 1/r = 0, which identifies the unpumped 
complex Maxwell-Bloch system as the purely nonlocal response, undamped limit (3?(7) = 
0, 1/r = 0) of the four- wave mixing model. 
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The second solution 1/r puts restrictions on the functions C and A. The condition 
on C, whose general solution is [H] 

2z/t = C + F(e~ 2t/T C), F arbitrary function, (33) 

restricts the expansion variable x t° only depend on the reduced variable £ = \/2ze~ t l T 
(the V2 is pure convenience) and therefore defines a reduction (z, t) — > £ of the PDE 
system to an ODE system written and studied in section [61 As to the restriction on A, 
which only makes sense for this £ reduction, it will be further examined in section [6j 

4.2 Conclusion of the test 

The result of the test provides the guidelines to be followed in order to obtain explicitly 
singlevalued solutions of the four-wave mixing model. These detailed guidelines, summa- 
rized in Table HJ are the following. 

• In the stationary case dt = 0, the test (not performed here) succeeds, therefore an 
eight-parameter singlevalued solution may exist. It has already been obtained in 
section [31 

• In the nonstationary, purely nonlocal response, undamped case (dt ^ 0, ^(7) = 0, 
1/r = 0), the system is equivalent to the unpumped complex Maxwell-Bloch system 
(I18p . integrable in the sense of the inverse spectral transform p], i.e. it admits N- 
soliton solutions, see section [5j 

• In the nonstationary, purely nonlocal response, damped case, no singlevalued so- 
lution exists unless the dependence on (z, t) is through the reduced variable £ = 
V2ze~ t/T . Then, a singlevalued solution may exist which depends on ten arbitrary 
parameters, we obtain it explicitly in section 16.11 

• In the nonstationary, arbitrary response case, whether damped or undamped, which 
includes the generic situation of the four-wave mixing, the structure of singularities 
is quite similar to that of the cubic complex Ginzburg-Landau equation (1121) (total 
differential order four, two irrational Fuchs indices, no movable logarithm [9]). Sin- 
glevalued solutions are locally represented by two Laurent series depending on eight 
(instead of ten as in the two previous cases) arbitrary functions, and the question 
to find closed form solutions a priori presents the same difficulty as for the CGL3 
equation. 

5 The unpumped Maxwell-Bloch system limit 

Since the pumped complex Maxwell-Bloch system (|18[) admits the Lax pair [TjJ] 

in which A is an arbitrary complex constant (the spectral parameter), the undamped four- 
wave mixing model with a purely nonlocal response then admits A^-soliton solutions, etc, 
as well as solutions in terms of the third Painleve function [251 12Q1 12] • 
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Table 1: Possible singlevalued solutions, according to time dependence, response (7) and 
damping (r). The reduced variable is £ = \f2ze~ l l T . 



d t M(j)1/t 


dependence 


solution 


Section 


= 






/(*) 


8-param 


m 




= 


= 


/CM) 


Maxwell-Bloch 


m 


7^0 


= 


7^0 


/(0 


10-param 


16.11 


7^0 


7^0 





/CM) 






7^0 


7^0 


7^0 


/CM) 







6 The dynamical case, reduction £ = (22:) 1 / 2 e *' 7 



The reduction (z,t) — > £ = (22) 1 / 2 e~*/ T (with an arbitrary origin for z and £) isolated 
by the Painleve test also exists for any value of 7 and we define it so as to preserve the 
definitions and (TT5T) . 



— 7^ 0, 7 arbitrary : < 
r 



f = e- 24/T -^/ m , r (0, Id(M) = e~ 2 * /T /d >r (e), 

^(z,t) = (l/2)e^-^(2^)- 1 /^ r (0, 
^•(z,t) = e-'/— '/%,(£), J = 1,4, 
{ A j (z,t) = e- t / T+i " t / 2 A j , I (0, j = 2,3. 



(35) 



It introduces one arbitrary real parameter uo. 

The intrinsic system (flBl - flHI) for I m ,Id,£ and the linear system for the amplitudes 



Aj reduce to 



1/ 

d£ m ' r 
d 



d_ * d 



%UJT 



d£ d£ d£ d£ 

[ iro = e 4t/r K(t) = / d 2 r + 4|J m , r | 2 . 

When compared to the travelling wave reduction (z,t) = z — ct,c ^ 0, 

I m (z,t) = e-^/ m . r (C), I d (z,t) = J d , r (C), £(M) = e-^ r (C), 
A,(^,t) = e-^/%. r (C), J = 1,4, 
A J (z,t) = e + ^ 2 A 3 ,(0, j = 2,3. 

d , . _ d . — ^ . , d „ 7 



(36) 



d cl d °d 

dC dC dC dC 

[ ir = K(t)=/ d 2 r + 4|J mir | 2 , 



-A (37) 
r c 



the two reduced systems (1561) and (1371) only differ by the evolution of the grating £ T . 
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6.1 Dynamical case, purely nonlocal response: general solution 

A direct computation of the conditions ( 1321) for both reduced ODE systems (l36i) and ( 1371) 
yields 

reduction (2z) 1/2 e' t/T : Q 2 = 0, Q 3 = cotC 3 , 
^ (7) = ° : ^ reduction z - at, c ? : Q 2 ee -L Q 3 = u (V* ~ \c 2 ' V) , (38) 

and the enforcement of Qj = makes both systems identical. Let us integrate the system 
(155]) with 3J(7) = 0, = 0. 

Thanks to the identity of the two systems (1361) and (1371) when the conditions = 
are enforced, the first integrals of the system (136|) for (I m ,r, Id,r, £ r ) can be generated 
systematically from the reduction X — cT of the Lax pair (1341) of the unpumped Maxwell- 
Bloch; this provides three first integrals, all real, 

K' = (7nlt) 2 Mdr + 4 l J m, r | 2 ) , 7 = «7nl, 7nl real, 



^1 = 7NLT (/m.r^r + ^m.r^rj , (39) 
3e = -7NLT-^d,r - |£r| 2 - 

Therefore 7d,r obeys a first order obtained by the elimination of £ v and / m ,D 

I' d / + 27N L r/ d 3 r - 12e / d 2 r - 2( 7NL r)- 1 ^/ dir + 4(K 2 + 3e ^)( 7NL r)- 2 = 0. (40) 

The general solution (J mr , Jd,r, ^r) of (|36ll) is singlevalued and expressible with the clas- 
sical functions p, (, a of Weierstrass, 

p' 2 = Ap 3 - g 2 p -g 3 = 4(p - e x )(p - e 2 )(p - e 3 ), p = C = (log a)'. (41) 

With the correspondence 

K' = g 2 - 12e 2 , K\ = -p\a)\ -2e = p(a), (42) 

the squared moduli and the gradient of their phases are doubly periodic functions, 

i r ,2 P' 2 (0 ~ P' 2 ( a ) r - p(a) |p|2 m , 

4(7nlt) 2 (-p(0 + p(a)) 7 NL r 

(arg/ m ,) --2^ _ , (arg£ r ) - - ^ — + 

e z(arg/ m , r - arg£ r ) = ^ - tp'(Q ^ 

27NLT|^r^m,r| ' 

the five constants of integration being eo, g2, #3 (actions), the origin of £ and the common 
origin of the phase of 7 mr and £ r (angles). 

The complex amplitudes themselves (I m;T ,£ T ,Aj) are also singlevalued functions and 
their expression, analogous to the complex amplitude of the traveling wave of the nonlinear 
Schrodinger equation, involves the a function of Weierstrass and is given in Appendix [HI 



1 When 3?(7) = 0,cj ^ 0, the ODE for I d:T has second order and is studied in [SJ Eq. (19.6)]. 
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An important particular case occurs for p'(a) = 0, all amplitudes then have constant 
phases. It proves convenient to first write this solution in complex form, in the symmetric 
notation of the Jacobi functions as introduced by Halphen [T6] . 



(x) = J p(x) - e a , a =1,2, 3, \imxh a (x) = +1, (44) 



L 



e 2 ^ha(0, £ 



-2iip 



e 2 ^h/3(e)h 7 (e), I m ,r 



i^ = — (-2h^(0 -3« 



7nlt. 



e = -- 



7nlt 



-2iip 



M0h 7 (£), 



« e 
i e 



+lf0 j 

-iipo 



012^(0+ fo i2h 7 (0) ; v4 2 ,r = ^oe" 
a 3 4h/3(0 + &34h 7 (£)), = «oe" 

(^12h/3(0+£l2h 7 (0), ^2,r = ^ e" 



¥0 I 



1<P0 I 



ai 2 h 7 (0+ &uM0)> (45) 
a 3 4h 7 (0 + hihp(£,)), 

-Al 2 h 7 (0 --Bijj h/3(0), 



3 , r = ^oe-^ (A 3 4h^(0 + £ 34 h 7 (0), ^ 4 < r = ^e + ^ (-A 3 4h 7 (0 - A* M0), 



7nlt 



with (a, /3, 7) an arbitrary permutation of (1,2,3) and the relations (!6lj) for the eight 
constants in Aj )r . 

In terms of the real Jacobi functions, the complex solution (1431) defines four bounded, 
physically admissible solutions (i.e. with positive square moduli for the amplitudes), in 
which the grating amplitude £ r is, respectively, a cn, dn, sd, nd function (with the usual 
notation k' 2 = 1 — A; 2 ), 



hi(C) = irk cn(r£, k), 112(6) = rk s n ( r 6, k), 113(6) = ir dn(r£, k), 
(a, /3, 7) = (1, 2, 3) : K' = r 4 , 6e = r 2 (l - 2k 2 ), 
(a, (3, 7) = (3, 2, 1) : K' = r A k\ 6e = r\k 2 - 2), 

hi(0 = irkk'sd(r^, k), 112(6) = r ^ cd(r£, fc), 113(6) = _ ' irfc'nd^, A;), 



(a,/5, 7 ) = (1,2,3): ^ = r 4 



6e c 



r 2 (l-2£; 2 ), 



(46) 



(47) 



{ (a, P, 7) = (3, 2, 1) : ^ = r 4 fc 4 , 6e = r 2 (fc 2 - 2). 



In these nine-parameter solutions, r, k are real, and A 12, A34 must be taken real in (l6Tj) to 
ensure that Aj and A, are complex conjugate. 

A second important case is the degeneracy from doubly periodic to simply periodic. 
The subcase p'(a) 7^ 0, which would correspond to the dark solitary wave 



le 



2iipo 



(fctanh(fc0 -m)e iK? , S r 



-2iipo 



(fctanh(fc6) +^)e" iK? , 



i 



7nlt 



■e aM (H; 2 + k 2 + ifc K tanh(fc6) - fc 2 tanh 2 (fc6))e^, 



'(/c 2 



ifcKtanh(fc6) - A; 2 tanh 2 (A;0)e" iK? , 



(48) 



7niT 



(-2A; 2 tanh 2 (fcO + 2fc 2 + /t 2 



is unphysical since the square modulus £ T £ T is negative. As to the subcase p'(a 
defines the bright solitary wave obtained from the long wave limit k 2 = 1 in (j46l 

hi (6) = h 3 (6) = irsech(r0, h 2 (6) = rtanh(r6), K = r \ K i = °, 6e = -r 5 



0, it 



(49) 
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with r, £ 0) V?o arbitrary. In this eight-parameter solution, r is real, and A12, A34, \x must be 
taken real to enforce the complex conjugation between Aj and Aj. 

Remark 1. Despite the similarity with the stationary value (|23|) for this bright profile 
of the grating amplitude, there is no limiting process yielding (|23|) from (|49|) . 

Remark 2. For those solutions displaying constant phases for the amplitudes, there 
must exist a value of the damped sine-Gordon variable u, Eq. ( TTUj) . able to represent the 
solution. Up to the numerous additive and multiplicative constants in (jTOl) and (145]) . this 
value is essentially given by 



e iu = h A (0, 

in which h.A and h a are related by the Landen transformation § 16.14.2] 



^loghA(0 



Mfl hc(Q 
iu(0 



ha(0, 



(50) 



(51) 



the correspondence between the elliptic invariants (e a , ep, e 7 ) and (e^e^ec) being de- 
tailed in [2], § 16.14.1]. For the trigonometric degeneracy (1491) . the value is 

e iu = rtanhr£, (52) 

and the Landen transformation reduces to the doubling of the argument with some shift, 



Vx : tanh x — 



tanhx 



— 2% sech 



2x + i- 



.7T 



tanh a; + 



tanhx 



2 tanh 



2x + i- 



.7T 



(53) 



7 Conclusion 

The four-wave mixing has been characterized by a lower- dimensional system of a deformed 
Maxwell-Bloch type. Then the three and only three possibly singlevalued limits of the 
four-wave mixing model have been determined and integrated. These consist of: (i) the 
stationary case for any r and 7; (ii) the limiting case 1/r = 0, 9^(7) = which is identified 
to the complex unpumped Maxwell-Bloch system; (iii) when ^(7) = 0, the reduction 
£ = \[2ze~ t l T to an ODE system. Those solutions which are localized (typically Jacobi 
bounded functions sn, cn, dn, cd, nd, sd [2, § 16.2]) should improve both the design of 
the physical devices to be manufactured and the confidence in the numerical simulations. 
As is often the case with methods based on singularities, the present study cannot rule 
out possible closed form but multivalued solutions. 

Moreover, the generic case 1/r ^ 0, §£(7) 7^ has been shown to display a structure 
of singularities, i.e. of possible closed form solutions, quite similar to that of the cubic 
complex Ginzburg-Landau equation. These solutions will be investigated in a forthcoming 
paper. 
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8 Appendix. Complex amplitudes of the integrable 
f reduction 

By elimination from ( 1361) . both fields S r and S r obey the same equation, 

^ - 2(p(f) - e )) V = 0, V = £■ (54) 

According to a classical result of Floquet, any linear differential equation with doubly 
periodic coefficients admits at least one solution which is doubly periodic of the second 
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kind [16J. The elementary unit of such doubly periodic functions of the second kind has 
been introduced by Hermite under the name element simple H(£, g, k) [161 v °l- II, p. 506], 



E(£,q,k) = ?§±A/ k -^, 



(55) 



chosen to have as only singularity a simple pole with residue 1 at the origin. Lame 
indeed proved that Eq. ([51]) admits the two solutions H(£, — a, 0) and H(£, +a, 0), which 
are generically linearly independent. Hence the complex amplitudes 



£ r = -ie 2i "H«,-a,0), /„,, 



27nl^ p{0 - p{a) 



-le 



-2iip 



H(£,a,0), I n 



27NLT 



-2itp 



P'(0 + P'(a) 

p{0 - pip) 



(56) 



H(£,a,0), 



in which the five constants of integration are eo, g 2 , #3, V^o and the origin of £. 

Given the values ( 1561) of £ r (£), £ r (£), eacn variable X, K of the linear system (flBi) also 
obeys a second order linear differential equation with doubly periodic coefficients, e.g., 



' d 2 1 p\j) - p'{a) d 

^e 2 2 P (o - P (a) de 



(P(0 - *>(<*)) U = 0, X = A 1>T . 



(57) 



This equation has the same features as (Jo4]) : unique singularity £ = of the Fuchsian 
type, Fuchs indices equal to —1,1, absence of logarithms in the general solution. A direct 
search for solutions of the elementary type ( 1551) provides the two solutions, generically 
linearly independent, 



X = H(f , +a/2 ± ft, 0), Y = H(£, -a/2 ± ft, 0), p{h) = p{a/2) - 2 



p /2 (a/2) 
p"(a/2) ' 



Taking account of the first integrals 

A\ >t A\ )T + A 2r A 2 ^ = constant, A 3r A 3r + A^A^ = constant, 
the general solution for the complex amplitudes can be parametrized as 



Ai, r = i ( a 12 H(£,+a/2 + M) + 612 H(£, +a/2 - ft, 0)) e +i ^°, 
A 2 ,r = io( a 12 H(£,-a/2 + ft,0) + b 12 H(f, -a/2 - ft, 0)) e - ^, 
A 3 , r = * ( a 43 H(e,+a/2 + a,0) + 6 43 H(£, +a/2 - ft, 0)) e+^°, 
^4,r = ^o( a 43 H(e,-a/2 + a,0) + 643 H(£, -a/2 - ft, 0)) e - *^, 
A ljr = * ( Aia H(e, -a/2 - ft, 0) + B 12 H(£, -a/2 + ft, 0)) e~^°, 
A 2 ,r = *o (-A12 H(£, +a/2 - ft, 0) - B 12 H(C, +a/2 + ft, 0)) e+^°, 
^3,r = *o ( A u H(£, -a/2 - ft, 0) + 5 34 H(£, -a/2 + ft, 0)) e~^°, 
A^ r = z (-^43 H(£, +a/2 - ft, 0) - 543 H(£, +«/2 + a, 0)) e+^°, 



(5? 



(59) 



(60) 
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with 



a±2 — cos cosh A12 e +iai2 , A i2 — cos /i cosh A 12 e" 
612 = i cos [l sinh A i2 e +i/312 , B 12 = i cos /i sinh A i2 e~ tl312 , 
a 34 = sin /i cosh A34 e +m34 , A 34 = sin /x cosh A 34 e~"* 34 , 
&34 = « sin/i sinh A 34 e +i/334 , _B 34 = i sin // sinh A 34 e~* 1 1 



-iai2 



(61) 



J(«12 " A2 - "34 + Am) = ±1) tan 2 ^ = ± 



sinh(2A 
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sinh(2A 34 ) 



The five additional integration constants are three of the four constant real phases a±2, 
P12, «34, P34, plus the two complex constants A12, A 34 . Finally, the conditions that Aj >r 
be the complex conjugate of Aj jT puts on A12, A 34 some constraints which depend on the 
choice for H(£, ±a/2 ±h,0), see text. 
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